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ABSTRACT 

The Chandra image of the merging, hot galaxy cluster 1E0657-56 reveals a bow shock propagating in front of 
a bullet-like gas cloud just exiting the disrupted cluster core. This is the first clear example of a shock front in a 
cluster. From the jumps in the gas density and temperature at the shock, the Mach number of the bullet-like cloud 
is 2-3. This corresponds to a velocity of 3000-4000 km s~' relative to the main cluster, which means that the 
cloud traversed the core just 0.1-0.2 Gyr ago. The 6-7 keV "bullet" appears to be a remnant of a dense cooling 
flow region once located at the center of a merging subcluster whose outer gas has been stripped by ram pressure. 
The bullet's shape indicates that it is near the final stage of being destroyed by ram pressure and gas dynamic 
instabilities, as the subcluster galaxies move well ahead of the cool gas. The unique simplicity of the shock front 
and bullet geometry in 1E0657-56 may allow a number of interesting future measurements. The cluster's average 
temperature is 14- 15 keV but shows large spatial variations. The hottest gas (T > 20 keV) lies in the region of 
the radio halo enhancement and extensive merging activity involving subclusters other than the bullet. 

Subject headings: Galaxies: clusters: individual (1E0657-56) — intergalactic medium — X-rays: galaxies 



1. INTRODUCTION 

Galaxy clusters form via mergers of smaller subunits. Such 
mergers dissipate a large fraction of the subclusters' vast kinetic 
energy through gas dynamic shocks, heating the intracluster gas 
and probably accelerating high energy particles (e.g., Sarazin 
2001 and references therein). Shocks contain information on 
the velocity and geometry of the merger. They also provide 
a unique laboratory for studying the intracluster plasma, in- 
cluding such processes as thermal conduction and electron-ion 
equilibration (e.g., Shafranov 1957; Takizawa 1999). Some ex- 
ploratory uses of X-ray data on cluster shocks were described 
by Markevitch, Sarazin, & Vikhlinin (1999). While many 
merging clusters exhibit recently heated gas (e.g., Henry & 
Briel 1995; Markevitch et al. 1999, Furusho et al. 2001 and 
references in those works; Neumann et al. 2001), so far only 
two candidate merger shock fronts were observed. One is a 
mild X-ray brightness edge, apparently a shock with a Mach 
number near 1, preceding the prominent "cold front" in A3667 
(Vikhlinin, Markevitch & Murray 2001). Another is a hot re- 
gion in front of the A665 core (Markevitch & Vikhlinin 2001, 
hereafter MV) which shows no clear density jump perhaps be- 
cause of an unfavorable viewing geometry. 

The Chandra observation of 1E0657-56 presents the first 
clear example of a cluster bow shock. This z = 0.296 cluster 
was discovered by Tucker, Tananbaum, & Remillard (1995) as 
an Einstein IPC extended source. From ASCA data. Tucker et 
al. (1998, hereafter T98) derived a temperature around 17 keV, 
making this system one of the hottest known (see also Yaqoob 
1999; Liang et al. 2000, hereafter LHBA). ROSAT data showed 
that 1E0657-56 is a merger (T98). It also hosts the most lumi- 
nous synchrotron radio halo (LHBA). 

Below we present results from the Chandra observation 
of 1E0657-56 performed in October 2000. We use Ho = 
lOO/i km s-'Mpc"' and f7o = 0.3, A = 0; 1' = 0.172/1-' Mpc 
at the cluster redshift. Confidence intervals are 90% for one- 

'AIsoatUCSD. 

^A combination of blank field observations normalized by the exposure time. 



parameter, unless specified otherwise. 



2. data ANALYSIS 

1E0657-56 was observed by ACIS-I at the focal plane tem- 
perature of -120°C for a useful exposure of 24.3 ks. To derive 
the gas temperature for a given region of the cluster image, the 
telescope and detector response were modeled as described in 
MV. The ACIS background rate did not vary during the expo- 
sure, but was higher than expected by a factor of about 1.3, 
most likely due to anomalous "space weather". This required 
special background modeling. To do this, we extracted a spec- 
trum from the ACIS-I region outside an r = 8.6' (1.5/;-' Mpc) 
circle centered on the cluster, which should be free of cluster 
emission. Point sources were excluded. The observed excess 
over the nominal background^ was well-modeled in the 0.7- 10 
keV band by the sum of two power laws with photon in- 
dices a = -0.6 and +3.0 (dominant below and above E 5 
keV, respectively) originating inside the detector, i.e., without 
applying mirror effective area and CCD efficiency to the model. 
Such a background anomaly in the ACIS-I chips is rare and not 
yet understood. We assumed that this component is distributed 
uniformly over the detector, and added it (normalized by solid 
angle) to the nominal background spectra. High-energy resid- 
uals in the overall cluster fit indicated that the normalization of 
the corrected background required an additional 10% increase 
(perhaps indicating some spatial nonuniformity of the excess), 
which we applied to the spectra from all cluster regions. 

To make a 0.5-5 keV image for the gas density analysis, 
we compared the observed background rate far from the clus- 
ter to the nominal model background (without the above addi- 
tional component) and derived a correction factor of 1.35 for 
this wide band, which was applied to the model background 
image. This is of course consistent with the above spectral cor- 
rection. A 10% background uncertainty was included in de- 
riving the confidence intervals for all quantities. This approxi- 
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Fig. 1. — {a) ACIS 0.5-5 keV image. Pixels are 4"; linear scale is for h = 0.5. {b) Grayscale shows optical ^-band image from ESO NTT 
(courtesy of E. Falco and M. Ramella). Contours, spaced by a factor of 2, show the smoothed ACIS image. The western part of the outer contour 
is approximately at the shock position. A number of well-matched point sources shows that the coordinates are accurate. 



mate background modeling is adequate for our present qualita- 
tive study. 

To minimize the effects of calibration uncertainties, the spec- 
tra were extracted in the 0.9-9.5 keV energy band, excluding 
the 1 .8 - 2.2 keV interval around the mirror Ir edge. 

3. RESULTS 

The ACIS image of the cluster is shown in Fig. la. It shows 
a "bullet" apparently just exiting the cluster core and moving 
westward. This subcluster was previously seen in the ROSAT 
data (e.g., T98), but the high resolution Chandra image makes 
its nature and direction of motion clear The bullet is preceded 
by an X-ray brightness edge that resembles a bow shock. To 
determine whether it is indeed a shock (and not a "cold front", 
e.g., Markevitch et al. 2000; Vikhlinin et al. 2001), below we 
deriv e the gas temperatures on both sides of the feature (see 
§3.1). Figure lb shows X-ray contours overlaid on an optical 
image. A subcluster of galaxies (e.g., Barrena et al. 2001) is 
seen leading the X-ray bullet, which is apparently swept back 
from the galaxies by the ram pressure of the ambient cluster 
gas. 

3.1. Temperature Map 

We first fit an overall cluster spectrum within r = 3' (0.5 /r' 
Mpc) as described in §^ using the Kaastra (1992) plasma 
model. We obtain T = 14.8;!;} j keV and an abundance of 
0. 1 1 ± 0. 1 1 solar, fixing the absorption at Nh = 4.6 x 10^° cm"^ 
as derived by LHBA from radio and ROSAT PSPC data. This 
temperature agrees with their ASCA + ROSAT fit of 14.5:1:^ 7 
keV and is consistent with the 17.4 ±2.5 keV ASCA fit by 
T98. UNh = 6.5 X 10^" cm'^ from Dickey & Lockman (1990) 
is used instead, we obtain T = 13.6 keV. When Nh is fit as 
a free parameter, it is consistent with both those values; the 
present low-energy ACIS-I calibration is not reliable for mea- 
suring Nh independently. Below we fix Nh = 4.6 x 10^° cm"^. 
Within the r = 0.5 Mpc aperture, the cluster's 0.5-5 keV 



(rest-frame) luminosity is 8.7 ±0.5 x 10'^* h ^ erg s ' and Lboi = 
2.3 ± 0.2 X lO'^^r^ergs-i. 

Figure 2 shows a temperature map made by dividing the 
cluster image into several regions and fitting their temperatures 
and abundances as above. Despite large uncertainties, the map 
shows that the gas outside the shock feature (region P) is cooler, 
or at least not hotter, than that inside (region S), confirming that 
it is indeed a shock front. The temperature at the tip of the bul- 
let is low (^ 7 keV) and is likely to have been the temperature 
of the subcluster. The hottest region of the cluster is its south- 
eastern X-ray brightness elongation. The optical image shows 
several large galaxies in that area, suggesting that this is the 
main merger site. As seen in Fig. 3, this also is where the radio 
halo is enhanced (LHBA). The halo also extends to the western 
shock front. A spatial correlation between the halo brightness 
and the local gas temperature (in addition to the general simi- 
larity to the X-ray brightness, e.g., LHBA; Govoni et al. 2001) 
was noticed by MV in two other merging clusters and supports 
the merger shock origin for the relativistic halo electrons (e.g., 
Tribble 1993). 

4. DISCUSSION 

4. 1 . The shock front 

The temperature map confirms that the western X-ray bright- 
ness edge is a shock front. We can derive its Mach number 
from either the temperature or density jump across the front us- 
ing the Rankine-Hugoniot shock adiabat. Figure 4a shows an 
X-ray brightness profile in a 120° sector centered on the bul- 
let's center of curvature and directed along its apparent mo- 
tion. There are two brightness edges whose shapes indicate 
spherical gas density discontinuities in projection. We fit this 
profile by the projection of a gas density model consisting of 
two power laws centered on the bullet representing the bul- 
let gas and the shock region, respectively. These are immersed 
in a /3-model centered on the main cluster (region j) represent- 
ing the outer, undisturbed gas. All components are spherically 
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symmetric around their respective centers; it is an adequate as- 
sumption in the sector of interest. Free parameters are the three 
slopes, two jump ampUtudes and two jump radii. A reasonable 
range of core radii and center positions for the /3-model was 
explored (obviously, these cannot be restricted by the fit) and 
found to have a small effect on our main interesting parameter, 
the density jump at the shock. Its confidence interval includes 
this modeling uncertainty. The observed temperature difference 
has a negligible effect on the derived density. 
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Fig. 2. — Projected gas temperature map (colors) overlaid on the 
X-ray brightness contours from Fig. lb. Temperatures for individual 
regions are shown in lower panel (region B is the tip of the bullet, 
unmarked in the map for clarity. This region appears smaller than the 
bullet in the contour plot because of smoothing). Errors in this figure 
are la. The horizontal line shows the cluster average temperature. 

The best-fit density model is shown in Fig. 4b; it describes 
the brightness profile in Fig. 4a well. The best-fit radial slopes 
are a w 0.15 for the bullet and a w -0.3 for the shock region; 
P « 0.7 assuming a core-radius of 125 kpc. The density 
jumps by factors of 3.8^[ q at the bullet edge and 3.2 ±0.8 at 
the shock front. Figure 4b also shows an approximate pressure 
profile (the density model times the temperatures from Fig. 2). 
The approximate pressure continuity at the first jump indicates 
that the bullet boundary is a "cold front", or contact disconti- 
nuity, similar to those recently discovered in other clusters. As 
expected, there is a large pressure increase at the shock front. 

A density jump of 3.2 ±0.8 at the shock corresponds to a 
Mach number M = 3.4 (> 2.1 at 90%) for a 7 = 5/3 gas and 
a one-dimensional shock (strictly speaking, the latter approxi- 



mation applies only along the direction of the motion, but we 
can use it as a qualitative estimate for our wide sector). The ob- 
served temperature jump from 8;!;3 keV to 18;!;5 corresponds 
toM = 2.1±l.l (again, for a qualitative estimate, we assume 
constant temperatures in regions P and S). These two indepen- 
dent derivations for the Mach number agree within their 90% 
uncertainties, and we conclude that M = 2-3. In the station- 
ary regime, the subcluster should move with the shock velocity. 
Such Mach numbers and the observed gas temperatures corre- 
spond to V ~ 3000-4000 km s"\ implying that the bullet has 
passed the center of the main cluster (^0.3li~^ Mpc away) just 
0.1-0.2 Gyrago. 

In principle, M could also be estimated from the Mach cone: 
the asymptotic angle if of the shock w.rt. the symmetry axis 
should satisfy sin (p = M"' . For M = 2-3,(p should be 20-30°, 
whereas the image suggests ip > 45°. A deviation of the ve- 
locity vector from the sky plane could widen the Mach cone in 
projection; however, a sufficiently large inclination angle also 
would smear the observed sharp brightness edges. An optical 
measurement of the subcluster relative line of sight velocity in 
the comoving frame, ^ 800 km s"' (Barrena et al.), combined 
with our value of M, also indicates only a small (10-15°) de- 
viation from the sky plane. The Mach cone relation assumes 
a uniform pre-shock medium and constant velocities, and the 
cone angle discrepancy is probably due to the cluster's radially 
declining density profile and deceleration of the bullet. 
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Fig. 3. — Radio halo brightness contours (reproduced from LHBA) 
overlaid on the temperature map from Fig. 2. The contours are at 
(3,6, 12, 18,24) xcr. 

4.2. The bullet subcluster 

The 0.5-5 keV luminosity of the bullet subcluster is ~ 5 x 
10'^^ erg s"', 10- 15% of the typical value for a cluster with 
the bullet's temperature; this fraction would change little if adi- 
abatic compression of the bullet by the surrounding hot gas is 
taken into account. The high gas density in the bullet (Fig. 4b) 
is of the order of that in the cooling flow clusters at compara- 
ble radii. This suggests that the bullet is a remnant of a density 
peak (cooling flow) that once was at the center of the merging 
subcluster, perhaps around its brightest galaxy (Fig. lb). The 
subcluster's less dense outer gas was probably shocked and ram 
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Fig. 4. — (a) ACIS 0.5-5 keV surface brightness profile in a 120° sector centered on the "bullet" and directed westward. The profile is extracted 
in elliptical segments parallel to the shock front; the r coordinate corresponds to an average radius within a segment. Errors in this figure are la. 
The histogram shows the best-fit model (see text). The corresponding unprojected model density profile along the symmetry axis is shown in panel 
(b), which also includes an approximate gas pressure model using temperatures in regions B, S, P from Fig. 2. Error bars on pressure correspond 
to errors in temperature. 



pressure stripped to the radius where the pressures are balanced 
(see Fig. Ab). Most of the stripped gas may reside in the north- 
south bar-like structure near the center of the main cluster seen 
in Fig. la (which may be a pancake in projection), since that is 
where the ram pressure on the moving subcluster was the high- 
est. Disruption of a cooling flow by a merger was considered 
theoretically by, e.g., Fabian & Daines (1991) and Gomez et al. 
(2000). 

The shuttlecock shape of the bullet clearly shows that it con- 
tinues to be actively destroyed by gas dynamic instabilities 
(e.g., Jones, Ryu, & Tregillis 1996). The gas being swept back 
from the cool bullet appears to be quite hot (region a in Fig. 
2) and clumpy (regions a, b), suggesting interesting physics at 
the interface of the two gases. This interface will be studied in 
detail using a longer Chandra observation planned for 2002. 

4.3. High overall temperature 

The existence of a few extremely hot clusters such as 
1E0657-56 have been used in the past to derive cosmolog- 
ical constraints under the assumption that the high tempera- 
ture indicates high virial mass. However, the results presented 
here show that the assumption of hydrostatic equilibrium in 
1E0657-56 can easily overestimate the mass by a significant 
factor due to the ongoing merger and a temporary increase in 
temperature (see, e.g., simulations by Ricker & Sarazin 2001 
and Ritchie & Thomas 2001). It is thus difficult to draw any 
strong cosmological conclusions without a better estimate of 
the cluster mass. 

4.4. Possible future measurements 

The unique shock and supersonic subcluster in 1E0657-56 
enable several interesting measurements. Since the gas density 
and temperature jumps at the shock are related, their accurate 
measurement may give the equation of state for the intraclus- 
ter plasma. Magnetic fields, a relativistic particle population 
with sufficient pressure, a significant lag between the electron 
and ion temperatures, and large nonthermal energy losses in a 



shock could all change the true or the apparent value of 7. 

From the velocity of the bullet and the density of the am- 
bient gas, we can derive the ram pressure on the bullet. The 
gas bullet appears to be pushed away from the subcluster's dark 
matter potential well just now, so the gravitational pull of the 
subcluster should equal the ram pressure. This may give an in- 
dependent estimate of the subcluster's total mass (Markevitch 
et al. 1999). 

The cluster also presents an interesting opportunity to con- 
strain the collisional nature of dark matter (e.g., Spergel & 
Steinhardt 2000; Furlanetto & Loeb 2001). There is a clear 
offset between the centroid of the bullet subcluster's galaxies 
and its gas. If one measures the location of the subcluster's 
dark matter density peak (e.g., from weak lensing or detailed 
modeling of the gas distribution), one may determine whether 
the dark matter is collisionless, as are the galaxies, or if it expe- 
riences an analog of ram pressure, as does the gas. 

5. SUMMARY 

The Chandra observation of 1E0657-56 presents a prototyp- 
ical example of a merger bow shock. The shock propagates in 
front of a cooler gas "bullet", apparently a remnant of a dense 
central region of the merging subcluster whose outer gas was 
stripped by ram pressure. The subcluster Mach number is 2-3 
and its velocity is 3000-4000 km s"'. Thus the subcluster tra- 
versed the main cluster core only 0.1-0.2 Gyr ago. The bullet 
is at the final stage of being destroyed by gas dynamic instabil- 
ities. Its gas lags behind the subcluster galaxies due to the ram 
pressure of the shocked gas. The hottest gas resides in a dif- 
ferent region of 1E0657-56 where additional merging activity 
occurs. The overall high temperature of 1E0657-56 is unlikely 
to represent its virial temperature due to the ongoing merger. 

We thank L. Grego, H. Tananbaum and the referee for useful 
comments. Support was provided by NASA contract NAS8- 
39073, grant NAG5-9217, the Smithsonian and the CfA fel- 
lowship program. 
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